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Milliarcsecond structure of water maser emission in two young 
high-mass stellar objects associated with methanol masers 
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ABSTRACT 

Context. The 22.2 GHz water masers are often associated with the 6.7 GHz methanol masers but owing to the different excitation 
conditions they likely probe independent spatial and kinematic regions around the powering young massive star. 
Aims. We compared the emission of these two maser species on milliarcsecond scales to deteimine in which structures the masers 
arise and to test a disc-outflow scenario where the methanol emission arises in a circumstellar disc while the water emission comes 
from an outflow. 

Methods. We obtained high-angular and spectral resolution 22.2 GHz water maser observations of the two sources G3 1 .58 1 +00.077 
and G33.641-00.228 using the EVN. 

Results. In both objects the water maser spots form complex and filamentary structures of sizes 18-160 AU. The emission towards the 
source G31. 581+00.077 comes from two distinct regions of which one is related to the methanol maser source of ring-like shape. In 
both targets the main axis of methanol distribution is orthogonal to the water maser distribution. Most of water masers appear to trace 
shocks on a working surface between an outflow/jet and a dense envelope. Some spots are possibly related to the disc-wind interface 
which is as close as 100-150 AU to the regions of methanol emission. 

Key words, stars: formation - ISM: molecules - masers ~ instrumentation: high angular resolution 



1. Introduction 

High-angular resolution observations of high-mass star forming 
sites in the 6.7 GHz methanol maser emission have revealed a di- 
versity of morphological structures from very simple composed 
of a few milliarcsecond (mas) wide spots to complex and ex- 
tended clouds of arcsecond sizes (Phillips et al. 119981 Walsh et 
al. [19981 M inier et al. "2000, Dodson et al. l2004l Bartkiewicz 
et al. l2009l Pandian et al. 201 1 ). Linear structures with a veloc- 
ity gradient are thought to be the signature of an edge-on disc 
or torus around young massive protostars (Norris et al. 119981 
Minier et al. l2000i l but there is evidence that some of them are 
associated with outflows (De Buizer l2003l l or even a shock prop- 
agating through a rotating dense molecular clump (Dodson et 
al.^004). Arc-like or ring-like morphologies of methanol emis- 
sion seem to prove a disc scenario (Bartkiewicz et al. 2009!, al- 
though most of those structures do not show signs of rotation 
but rather inflow or/and outflow dominate (van Langevelde et 
al.|2pi0j, it is then postulated that methanol masers arise in the 
shock interface between the large scale accretion and a circum- 
stellar disc. 

While the 6.7 GHz methanol maser line is radiatively 
pumped in warm (T ~ 150 K) and dense (n < 10** cm"^) re- 
gions (Cragg et al. fmUS^. Sobolev et al. ^19971, the 22 GHz wa- 
ter maser emission is collisionally pumped and probes denser 
(n > 10^ cm-^) and hotter (r ~ 400 K) gas (Elitzur et al. [T989b . 
The water maser is often excited in strong shocks, driven by 
young low-mass and high-mass (proto)stellar objects, on a work- 
ing surface between a outflow and a dense envelope. A number 



of morphologies of this emission is observed from collimated jet, 
wide-angle flows, expanding shells to equatorial flows (Goddi 
et al. 2007; Moscadelli et aL l2007l IIOTTI Sanna et al. I2010al 
2010b i Torrelles et al. 201 l l. One or more types of those maser 
structures occur simultaneously in star forming regions on scale 
sizes of a few mas to several arcsec and is possibly related to the 
geometry of the envelopes or the mechanisms driving the out- 
flows. 

Our recent VLA survey of 22 GHz water masers in a sample 
of 3 1 methanol sources has yielded 22 detections (Bartkiewicz et 
al. l201 11 1 where both maser species are excited by the same cen- 
tral objects. We noted that a distinct group of methanol sources 
with ring-like structure show either no associated water masers 
at all or water masers that are distributed orthogonally to the ma- 
jor axis of the ring. It is argued that the methanol maser structure 
traces a circumstellar disc/torus around a high-mass young stel- 
lar object while the water masers originate in outflows. As the 
VLA data are of moderate angular and spectral resolutions of 
~1" and 0.65 kms ', respectively, we undertook the VLBI ob- 
servations to examine a disc-outflow scenario in the two bright- 
est water maser targets, which show the arc-like or ring-like mor- 
phologies that are characteristic of methanol masers. 

The source G3 1.581 +00. 07'fl has been recognized as a mas- 
sive young stellar object based on 6.7 GHz methanol maser ob- 
servations (Szymczak et al. 2000). The methanol maser spots are 
distributed along an arc or ring of 217 mas size (Bartkiewicz 



' The names of two targets follow the Galactic coordinates of the 
brightest methanol maser spots derived by Bartkiewicz et al. ( i20098 . 
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et al. l2009l l. which corresponds to ~1200AU for the assumed 
near kinematic distance of 5.5 kpc (Reid et al. l2009l l. The 
detection of an infrared source of bolometric luminosity of 
3xlO'*Lo (Urquh art et a l. l20Tll l. 22 GHz H2O maser emission 
(Bartkiewicz et al. l2011l l. 1665 and 1667 MHz weak OH masers 
(Szymczak & Gerard I2004I I. millimeter molecular thermal lines 
(Szymczak et al.'200T; Urquhart et al. 2008 1, and 5 and 8.4 GHz 
continuum emission at location of 31.582+00.075 (i.e., 9" apart 
from the maser arc/ring) (Urquhart et al. 2009 , Bartkiewicz et 
al. l2009l) indicates that this is a cluster of recent star formation. 

The source G33.641 -00.228 detected in the 6.7 GHz 
methanol line (Szymczak et al. 2000) has an arc spot distri- 
bution of length 630 AU (Bartkiewicz et al. |2009l l for the as - 
sumed near kinematic distance of 3.8 kpc (Reid et al. l2009l l. 
which seems to be more likely than the far kinematic distance 
(Zhang priv.comm.). The site also contains water and OH masers 
(Bartkiewicz et al. 1201 11 Szymczak & Gerai'd l2004t but no 
8.4 GHz continuum emission was detected with an upper limit 
of 0.15Jy (Bartkiewicz et al . 120091 1. The detection of radio re- 
combination lines at about 103 km s ' (Anderson et al.'201 1) and 
a multi-feature '•'CO line spectrum (Urquhart et al. 2008i illus- 
trates the complexity of this molecular cloud and its clustered 
star formation. 



2. Observations and data reduction 

The EVr^ observations of G3 1.58 1+00.077 and 
G33. 641 -00.228 with the antennas at Jodrell Bank, Effelsberg, 
Medicina, Metsahovi, Onsala, and Yebes, were carried out at 
22.23508 GHz on 2010 October 30 for 8 h (the project EB047). 
The tracking phase centres were estimated from the VLA 
survey for the water maser spots that were located nearest to 
the methanol emission in each target (Bartkiewicz et al.2011) 
at a=18''48™4P.951, (5=-01°10'02'.'578 and a=18*'53™32^563, 
5=+00°31'39'.'130 (J2000), respectively. A phase-referencing 
scheme was applied with a reference source J 185 1+0035 (from 
the VLA calibrator catalogue), using a cycle time between the 
maser and phase-calibrator of 60 s+90 s. This strategy yielded 
2 h on-source times for each target. The bandwidth was set to 
8 MHz and data were correlated with the Mk IV Data Processor 
operated by JIVE with 1024 spectral channels. The resulting 
spectral resolution was 0.1 km s"'. The velocity was measured 
with respect to the local standard of rest (LSR). 

The data calibration and reduction were carried out with 
NRAO's Astronomical Image Processing System (AIPS) using 
standard procedures for spectral line observations. We used the 
Effelsberg antenna as a reference. The amplitude was calibrated 
by performing measurements of the system temperature at each 
telescope and applying the antenna gain curves. The parallac- 
tic angle corrections were subsequently added to the data. The 
source 3C454.3 was used as a delay, rate, and bandpass cali- 
brator. The phase-calibrator was imaged and a flux density of 
190mJy was obtained. The maser data were corrected for all 
Doppler effects and self-calibrated using the brightest and most 
compact maser spot. Finally, naturally weighted maps of spectral 
channels were created in the velocity range where the emission 
was seen in the scalar-averaged spectrum. For imaging, the re- 
sulting synthesized beam was 1 .0x2.4 mas^ with a position angle 
(PA) of -38°, while the pixel separation was 0.2 mas in both co- 
ordinates. The rms noise levels (lo-^ms) in line-free channels was 

^ The European VLBI Network is a joint facility of European, 
Chinese, South African and other radio astronomy institutes funded by 
their national research councils. 



typically lOmJybeam The weakest detected maser spot had 
a brightness of 94mJy beam"' that is more than 9o-nns- 

The positions of water maser spots in all channel maps 
were determined by fitting two-dimensional Gaussian models. 
The formal fitting errors resulting from the beamsize/signal-to- 
noise ratio were smaller than 0.1 mas. To determine the posi- 
tion accuracy of registered maser spots, we need to consider 
the following factors (Diamond et al. 2003 ): i) the uncertainty 
in the phase-reference source position of 2maS ii) the an- 
tenna positions, where the claimed accuracy of ~ 1 cm corre- 
sponds to 1 mas in RA and 2 mas in Dec; iii) the phase trans- 
fer over the separation between targets and the phase-calibrator: 
l°.92 (between G3 1 .58 1 +00.077 and J 1 85 1 +0035) and 0^45 (be- 
tween G33.641 -00.228 and J1851+0035). These caused poten- 
tial phase-solution transfer errors corresponding to 1 mas in RA 
and 2 mas in Dec for G3 1 .58 1 +00.077 and to 0.2 mas in RA and 
0.4 mas in Dec for G33. 641 -00.228, respectively. In total, the 
absolute position accuracy (IcTpos) is 2.5 mas in RA and 3.5 mas 
in Dec for both targets. 

3. Results 

Water maser emission towards both targets was mapped, after 
successful phase-referencing, in the areas of 5"x5" and the 
velocity range explored by Bartkiewicz et al. (20Hi using the 
VLA. The EVN maps of water maser spots (Figs. [T]and|3]l are 
overlayed on the VLA maps and the 6.7 GHz methanol maser 
distributions obtained with the EVN (Bartkiewicz et al. 2009}. 
The Spitzer GLIMPSE map^ of the 4.5 yum-3.6 /im emis- 
sion excess are added. Following the procedure described in 
Bartkiewicz et al. (2009), the water maser spots are analysed 
to identify individual velocity-coherent maser clouds. Their ba- 
sic parameters, such as positions, ARA, ADec, LSR velocities, 
Vlsr, and intensities of the brightest spot, Sp, of each cloud 
are listed in Tables [T] and |2] The water and methanol maser 
spectra of single clouds are combined in Figures |2] and with 
overlays of individual Gaussians, if emission was seen in at 
least three consecutive channels. A Gaussian analysis of indi- 
vidual clouds is carried out and the fitted flux amplitude, Sfit, 
the full width at half-maximum, FWHM, the projected extent 
between the most separated single spot centres within a cloud, 
Lproj, and the velocity gradient, Vgrad are listed in Table |3] The 
lower limit to the brightness temperature, Tb, of each cloud 
is also calculated according to Eq. 9-27 of Wrobel & Walker 
(I1999I I. For comparison purposes, the same analysis is done for 
the 6.7 GHz methanol maser data (Table |3]l obtained in 2007 
June (G3 1.58 1+00.077) or 2003 June (G33.641 -00.228) with 
the EVN. Below we comment in more detail on each source. 

3.1. G31. 581 +00.077 

A total of 91 water maser spots were detected. They are con- 
centrated in two distinct regions (Fig. [T}. The south-east (SE) 
region of size 50x30mas^ containing weak (<1.35 Jy) emission 
in the velocity ranges from 90.1 to 93.9 km s"' and from 100.0 
to 103.3 km s"' is located close to the phase centre. The second 
region (NW) of middle velocity emission (96. 8-100. 7kms"') 
and flux density of 2.6-16.5 Jy is located offset by ~5" to the 
north-west (PA= -60°) of the phase centre. 

The SE region is divided into the cluster of clouds 7, 2, and 
3 separated by 50 mas from the cluster of clouds 6, 7, and 8 

3 The GSFC ICRF2 VLBI Source Position Catalog. 
http://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/ 
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(Fig. [TJ. The clouds in the SE region show Gaussian velocity 
profiles with FWHM linewidths of 0.54-0.94kms Their pro- 
jected sizes are 0.07-0.47 mas, which correspond to 0.9-2.6 AU 
for the adopted distance of 5.5 kpc. The brightness temperature 
is 0.24-1.38x10''K. We note that clouds 3 and 6 show the ve- 
locity gradient of 1.3-1.8kms"'mas"' oriented along PAs of 1° 
and -23°, respectively. It is remarkable that the two clusters of 
clouds (1, 2, 3 and 6, 7, 8) appear very close in space (Table [1} 
forming arc-like filaments of projected length of 4 and 3.2 mas 
corresponding to 22 and 18 AU, respectively, and elongated at 
PA=-58°. The closest clouds 3 and 6 are separated by 285 AU 
and a velocity difference of 6.7kms The mean linear sep- 
aration between the two clusters along PA=65° is 310±8 AU, 
while the velocity spread is 13.2kms The velocities of the 
two filaments are approximately symmetric with the regard of 
the systemic velocity of 96.0 km s"' (Szymczak et al. 2007). The 
two SE filaments are likely signatures of flattened shock sur- 
faces (Torrelles et al. ,200r ). We found that the water cloud 7 
coincides within +0.1 km s ' with the methanol maser cloud k, 
while the spatial separation is 24.4 mas which corresponds to the 
projected distance of 134 AU. Although these two lines were ob- 
served within a time span of three years their separation seems to 
be real. This confirms that the water and methanol masers probe 
different parts of the environment of young massive stars (e.g., 
Beuther et al. l2002l Sanna et al. l2010al |2010b") because of the 
different pumping mechanisms affecting both species (Elitzur et 
al.[l989;Cragg et al. l2005T l. 

The NW region of elongation of 26 mas is composed of two 
clouds 4 and 5 of brightness temperatures of 2.7-17.1xlO^K 
(Table O. Their velocity profiles are obviously asymmetric and 
nicely fitted by the sums of three and two Gaussian components 
(Fig.|2]l of FWHM Hnewidth of 0.5-0.9 km s'. The velocities of 
the NW water maser emission largely coincide with those of the 
methanol masers but their linear separation is about 27000 AU. 
The map of 4.5 //m- 3.6 fim emission excess (Fig. [TJ, which 
possibly traces shocked molecular gas in outflows from massive 
stars (e.g., Davies et al. l2007l l. shows that the SE and NW maser 
clusters lie in a very complex region. Both maser clusters are 
likely associated with different powering sources. 

The comparison of our maps with those obtained about 
14 months earlier with the VLA in a CnB configuration 
(Bartkiewicz et al. 1201 11 1 implies that the cloud 1 coincides to 
within 13 mas, the cloud 4 to within 56 mas, and the cloud 7 to 
within 64 mas with the matching VLA spots. The emission seen 
with the VLA from three components located ~l-2" eastward 
of the clouds 4 and 5 is not detected. We also note that the water 
maser components detected have very similar amplitudes in both 
EVN and VLA observations. 

3.2. G33.641 -00.228 

Twenty-one water maser spots were detected in the velocity 
range of 54.4-57.6 km s"' over the north-south elongated area 
50x280 mas^ close to the phase centre (Fig.|3]l. Using the above- 
mentioned procedure, we found four maser clouds where the 
emission is seen in at least three contiguous spectral channels. 
There are also two other clouds (5 and 6) that do not obey this 
criterion but seem to be real (Table |2}. All but one of the cloud 
are in the southward cluster at a distance ~270 mas from the 
phase centre. This filament cluster of size 41.6 mas is aligned 
along a PA of 78°. For the adopted distance of 3.8 kpc, the cor- 
responding linear scales are 1020 AU and 157 AU. The clouds 
1-4 in the southern cluster have Gaussian profiles with FWHM 
linewidths of 0.25-0.88 km s"'. We note that none of the water 



Table 1. List of 22.2 GHz water maser clouds as found in 
EVN observations towards G31.581-hOO.077. The (0,0) point 
corresponds to the position of the brightest methanol maser 
in this source (Bartkiewicz et al. I2009] l: RA=18H8'M1';94108, 
Dec=-01°10'02'.'5281 (J2000). 



Cloud 


Vlsr 


ARA 


ADec 


Sp 


Number 




(kms-') 


(mas) 


(mas) 


(mjybeam ') 


of spots 


1 


102.9 


137.78 


-43.02 


1008 


10 


2 


101.7 


136.79 


-42.58 


233 


4 


3 


100.2 


134.57 


-40.69 


472 


7 


4 


99.7 


-4008.86 


2763.28 


16534 


20 


5 


97.5 


-4003.61 


2789.88 


2577 


18 


6 


93.5 


88.23 


-64.04 


363 


7 


7 


91.7 


86.79 


-63.97 


1000 


14 


8 


90.7 


89.65 


-65.48 


1333 


11 



Table 2. List of 22.2 GHz water maser clouds as observed us- 
ing EVN towards G33.641 -00.228. The (0,0) point correspond 
to the brightest methanol maser in this source (Bartkiewicz et 
al. |2009] l: RA=18''53'"3Z563, Dec=+00°31'39'.'180 (J2000). 



Cloud 


Vlsr 


ARA 


ADec 


Sp 


Number 




(km s"' 


(mas) 


(mas) 


(mjybeam 


of spots 


1 


58.1 


87.77 


-271.39 


174 


3 


2 


57.3 


17.79 


-8.96 


606 


7 


3 


56.1 


46.98 


-279.33 


157 


4 


4 


54.8 


48.62 


-277.37 


263 


3 


5 


55.1 


48.65 


-277.49 


118 


2 


6 


54.5 


48.69 


-277.40 


217 


2 



clouds have an internal velocity gradient. The linear size of the 
individual clouds is 0.3-l.OAU and the brightness temperature 
is 0.12-0.27xlO^K (Tabled. 

The strongest emission (Tb = 0.63 x lO^K) comes from cloud 
2, which is a blend of two Gaussian components (Fig.|4]i. Cloud 
2 is ~20 mas (75 AU) away from the brightest methanol maser 
clouds a and r and differ in terms of velocity by 1 .5 km s"' . This 
is probably the first such tight association of methanol and water 
masers. We note that the emission of cloud 2 is blue-shifted by 
4.2 km s"' and that of cloud a is red-shifted by 1 .2 km s"' relative 
to the systemic velocity of 61.5 km s"' (Szymczak et al. 2007]l. 
It is therefore possible that these clouds signpost a shock front 
in which the water emission originates behind the front, while 
the methanol emission (cloud a) appears outside of the shock 
interface in the infalling gas. 

The image of the 4.5 //m-3.6 jjm emission excess (Fig. 13) 
shows that the water maser emission is located 1'.'2 to the south- 
east of single mid-infrared object. This is likely the source pow- 
ering the outflow along a PA of 165° traced by the water masers. 

The positions of the clouds 2, 3, 4 and 6 differ by 
about 92 mas from the VLA positions of the matching clouds 
(Bartkiewicz et al. l201 11 1, which is well within the ~150mas ac- 
curacy of the CnB configuration VLA data. A weak (<0.9Jy) 
red-shifted emission in the velocity range of 83.8-85.1 km s"' 
seen with the VLA towards cloud 2 (Fig. [TJ was not recovered 
in the VLBI observation. The brightness of the emission from 
cloud 2 determined in the EVN observation is only slightly lower 
than that measured with the VLA. In contrast, the maser spots 
from the southern region (clouds 1, 3-6) are about ten times 
weaker in the EVN maps that may suggest the presence of dif- 
fuse and/or highly variable emission. 
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G31.581+00.077 



o) 




b) 



4000 
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-4010 




Fig. 1. G31.581-h00.077. a) Spitzer GLIMPSE 4.5 yum-3.6 /jm excess image overlaid with the water maser positions (yellow aster- 
isks) from the EVN observation and H II region (yellow cross) detected by Bartkiewicz et al. (2009). b) Shows the distribution of 
water masers observed with the EVN (filled squares) from this paper and the VLA (open squares) from Bartkiewicz et al. (201 1}. 
The colours of the symbols relate to the LSR velocities as indicated on the right-hand side of the plot. The origin of the map is the 
position of the brightest 6.7 GHz methanol maser spot (Bartkiewicz et al . 120091 1 (see also Table[T]i. c) Shows an enlargement of the 
north-western (NW) water masers. d) Shows an enlargement of the south-eastern (SE) water masers together with the distribution 
of the 6.7 GHz methanol masers, marked by open circles, from Bartkiewicz et al. (2009). In (c) and (d), the arrows represent the 
velocity gradients (from blue- to red-shifted LSR velocities) detected within individual cloud. The black numbers and letters corre- 
spond to the clouds with internal velocity gradients, while the grey ones correspond to the clouds without internal velocity gradients 
(Table|3]l. 



4. Discussion 

4.1. Maser environments 

The Spitzer GLIMPSE map (Fig. [1]) clearly indicates that the 
source G3 1.58 1+00.077 lies in a large (35"x25") complex area 
with an 4.5 pm emission excess in the form of clumped and dif- 
fuse structures. This extended emission is proposed to be a tracer 
of shocked gas from regions where outflowing gas interacts 
with the surrounding medium (Davis et al.'2007 , Cyganowski et 
al.[2009). The '^CO(1-0) spectrum obtained with a 46" beam 
(Urquhart et al. l2008] l towards the position 31.5808+0.0757 has 
two Gaussian components centred at 96.2 and 109.8kms A 
weak emission line is clearly detected in the range of 85 to 
105 km s"', i.e., in the wings of the first component suggests 
the presence of outflows. A signature of ordered motions also 
appears in the HCO^(l-O) spectrum (Szymczak et al. l2007l l. 
where in this slightly asymmetric profile the strongest emission 



is blueward of the source systemic velocity and can be inter- 
preted as the result of infall (e.g.. Fuller et al. 2005 1. We then 
find 9'.'2 east of the SE water maser an H II region of flux den- 
sity 10.3 and 15m Jy at5 GHz (White et al. l2UU5]) and 8.4GHz 
(Bartkiewicz et al.'2009), respectively. This continuum source 
coincides to within 0'.'6 with a clump of 3.6 jum-4.5 pm emis- 
sion excess (Fig. [T]). The NW water masers are located at the 
borders of two clumps of excited gas, and their powering source 
is unclear although it seems unlikely that the SE and NW masers 
are associated with the same central object. A comparison of 
the present EVN maps with those obtained with the VLA 14 
months earlier (Bartkiewicz et al. l201 1) implies that there is a 
lack of water masers eastward of the NW clouds (Fig. [1). It is 
likely that these are diffuse and low intensity masers resolved 
out with the EVN beam. The spectrum obtained with a 40" 
beam sometime between November 2009 and December 2010 
(Urquhart et al. l201 11 1 indicates that the peak flux density of 
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4000 



3000 



G31.581+00.077 

22.2 GHz woter masers 

6.7 GHz methanol masers 



-£ 1000 




90 92 



96 98 
Vlsr (km s-') 



104 98 99 100 101 



Fig. 2. Individual component spectra with Gaussian velocity profiles of 22.2 GHz water masers (red squares and lines) and 6.7 GHz 
methanol masers (blue circles and lines) towards G3 1.58 1+00. 077 detected using the EVN in 2010 and 2007, respectively. The 
numbers and letters correspond to the cloud labels as given in Table[T]and Fig.[T] The grey lines show individual Gaussian profiles 
fitted to the blended features. The clouds with only two single spots are not marked to improve the clarity of the figure (i.e., clouds 
b, e, g). 



G33. 641 -00.228 



°) 



b) 
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Fig. 3. Same as Fig.[I]but for G33.641 -00.228. 
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-278 
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ARA (mas) 
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125 Jy at 99.7 km s"' is about one order of magnitude higher than 
that measured with the EVN. We suggest that the water masers 
of G31.581-hOO.077 are produced by different young stellar ob- 
jects in a complex region composed of at least a few high-mass 
stars well-signposted by the radio continuum and the methanol 
and water masers. 

The water maser in G33. 641 -00.228 lies only I'.'l south- 
eastward of the maximum of the 3.6 //m-4.5 fim emission excess 
of cometary-like morphology (Fig. |3). The HCO^(l-O) spec- 



trum obtained towards the position 33.648-0.224 (Szymczak et 
al. l2007l l shows a small dip near 61.5kms ' and a slight red- 
ward asymmetry that is indicative of outflow motions. The less 
optically thin H'^CO^(l-O) transition is detected as marginal 
red-shifted emission, which is consistent with an outflow. In 
the region where the OH 1665, 1667, and 1720 MHz masers 
were detected (Szymczak & Gerard l2004) and all of them peak 
at 60.2 km s"'. The strongest 1720 MHz emission and a broad 
(7.4kms ') 1667MHz absorption profile near 56.2kms ' are 
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54 55 56 57 58 59 60 61 62 63 64 
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Fig. 4. Same as Fig.|2]but for G33. 641-00. 228. The scale on the left side describes the water maser intensity, while the right side 
corresponds to the methanol maser intensity. The clouds with only two single spots or without Gaussian characteristic are not 
marked for the clarity of the figure (i.e., clouds 5, 6, e, g, i, n, p). 



indicative of shock fronts and a continuum background source, 
respectively. No continuum emission at 8.4 GHz was detected 
with a 3cr„„s limit of 0.15 mJy beam"' (Bartkiewicz et al. 2009). 
A weak (~20mJy) 8.6 GHz hydrogen radio recombination- 
line near 102. 9 km s" ' detected at position G033.645-0.227 
(Anderson et al. l201 11 1 differs so greatly in terms of velocity from 
the water masers that this may be a chance projection that is 
unassociated with the methanol and water masers. We argue that 
the water masers of G33. 641 -00.228 as well as the associated 
methanol and hydroxyl masers are excited by individual high- 
mass stars. 

4.2. Physical parameters of circumstellar medium 

Lower limits to the brightness temperature of the SE water maser 
components associated with G3 1.58 1+00. 077 are always lower 
than 1.4xlO'^K, while the measured linewidths range from 0.54 
to 0.94 km s ' are quite commonly in the star-forming sources 
(e.g., Goldreich & Kwan il974| Surcis et al. l201 lall20Tft] l. Since 
the kinetic temperature of the masing gas might be expected 
to be about 400 K (Elitzur et al. l 19891 1. the intrinsic thermal 
linewidth given by Avfwhm = 2.35482 x yfkTiJm, where Tk 
is a kinetic temperature, k the Boltzmann constant, and m is a 
molecular mass should be ~1 kms"'. This value is larger than 
the observed linewidths and suggests that the masers are unsatu- 
rated. In the same volume of gas where the water maser cloud 7 
and the methanol maser cloud k appear to coincide (Fig. [TJ, the 
methanol thermal linewidth would be 0.75 km s ', whereas the 
measured value is 0.30km s '. The narrowing of the line pro- 
file is expected when the maser is unsaturated. Detailed calcula- 
tions have shown that in one of the strongest known methanol 
sources NGC7538 about 92% of the components are unsatu- 
rated (Surcis et al. 120 lib) . In the source G33.641 -00.228, the 
observed linewidths of methanol components are generally nar- 
rower than in G3 1 .58 1 -1-00.077, which suggests that the maser is 
also unsaturated in this region. 

The mas-scale velocity gradients in both lines are ob- 
served in G3 1.58 1+00. 077. We note that the velocity gradi- 
ent at 22 GHz is about one order of magnitude higher than at 
6.7 GHz (TableO. Extensive discussion of the velocity gradients 
for methanol masers by Moscadelli et al. (1201 lb) suggests that 
there is a kinematical interpretation of their origin. We note that 
the methanol maser gradients in both targets can reflect the or- 
dered motions on scales of 10-60 AU. Large velocity gradients 
for water masers in G3 1.58 1+00. 077 suggest that they are re- 



Table 4. Parameters derived by fitting the kinematics of the ro- 
tating and expanding disc model. The signs + and - of the ro- 
tation and expansion velocities refer to the clockwise or anti- 
clockwise rotation and outflow or inflow for positive i. Both ro- 
tation and flow are reversed in the case of negative i. For each 
source, both signs together with the sign of i could be reversed 
since our model does not give the directions unambiguously. In 
the case of G33. 641 -00.228, two model fits are presented as de- 
scribed in Sect. 4.3. 



Source 


(kms-') 


(kms-') 


V' 

sys 

(km s-') 


i 

n 


xl 


G31.581+00.077 














0.90 


1.71 


98.54 


77.5 


134 


G33.641-00.228 












Dl 


0.84 


-1.05 


61.00 


30.0 


184 


D2 


3.72 


1.74 


58.40 


-30.0 


222 



Systemic velocity. 



Table 5. Parameters derived from fitting the biconical out- 
flow model of Moscadelli et al. (2000). In the case of 
G3 1.58 1+00. 077, two model fits are presented as described in 
Sect. 4.3. 



Source 


Xo 






PA 


e 








(mas) 


(mas) 


(kms-') 


D 


D 


n 




G31.581+00.077 
















Ol 


72 


68 


-77.0 


-16 


15 


83 


1.7 


02 


110 


-46 


7.7 


67 


17 


45 


0.7 


G33.641-00.228 


15 


5 


15.2 


16 


9 


67 


0.1 



' Constant velocity of a maser spot, + for a direction radially outward 
from the central star. 



lated to the outflow motions with velocities larger than 50 km s 
(see Sect. 4.3). 
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4.3. Kinematic models 

The two sources investigated in this study possibly belong to a 
group of objects where a ring-like or arc-like methanol maser 
distribution traces a circumstellar disc/torus around a high-mass 
young stellar object, whereas the water maser distribution is or- 
thogonal to the major axis of the methanol structure (Bartkiewicz 
et al.20ir). Furthermore, these objects are not associated with 
detectable continuum emission at cm wavelengths (Bartkiewicz 
et al. 2009) and may represent an early stage of evolution. The 
new EVN observations of both sources have shown that all of 
the individual water maser spots detected previously with the 
VLA (Bartkiewicz et al. l201 11 1 unfold into complex and filament 
structures of sizes 18-160AU. To examine a disc-outflow sce- 
nario in the two sources, we used a model of a rotating and ex- 
panding thin disc (Uscanga et al. l2008l l for the methanol masers 
and a model of the outflow (Moscadelli et al.'2000) for the wa- 
ter masers. Detailed descriptions of the modelling are given 
in the above cited works, and elsewhere for methanol masers 
in discs (Bartkiewicz et al. |2009l and outflows (Bartkiewicz et 
al. l2006l l. We note that in the following we analyse the new wa- 
ter maser data with the methanol maser data obtained 3-7 years 
earlier (Bartkiewicz et al. 2009 1. However, our inspection of a 
few sources published so far and systematic EVN monitoring 
demonstrate that the overall methanol maser morphologies are 
stable on a scale time of 6-8 years (Bartkiewicz et al., in prep.). 

In Tables 4 and 5, we summarize the best-fit values of both 
models. For the rotating and expanding thin disc, the rotation 
(Vrot), expansion (Vexp), and systemic (Vsys) velocities as well as 
the inclination angle, i, i.e., the angle between the line-of-sight 
and the normal to the ring plane defined to be i= acos(^) is given 
for each source. The solutions were based on the minimization 
of the function given by eq. 8 in Uscanga et al. (.2008 J . The 
outflow model for water masers is characterized by the vertex of 
the cone (Xq, Yq), the x-axis coinciding with the projection of 
the outflow on the plane of the sky at the position angle PA, the 
inclination angle between the outflow axis and the line-of-sight 
(i.e., the z-axis) and the opening angle of the outflow/cone 20. 
The systemic LSR velocities, Vc, were assumed to be the same 
as given in Sect. 3. The;t'^ function was assumed to be expressed 
as in Eq. 3 of MoscadelH et al. (2000) . 

The methanol emission in G3 1.58 1-1-00.077 is most accu- 
rately reproduced by the model where the powering source is 
at the centre of the best-fit ellipse and the rotation velocity of 

0. 9 km s"' is a factor of two lower than the expansion veloc- 
ity (Tab. m Fig. |5j. These values are typical of the class of 
ring-like methanol masers (Bartkiewicz et al. 12009) and sug- 
gest that the methanol structure of ~1000AU diameter cannot 
be interpreted as a Keplerian disc. It is instead proposed that 
the methanol masers arise on the interface region between a 
large-scale accretion flow and a stellar disc (van Langevelde et 
al. l2010l Torstensson et al. l201 11 ). We also note that the best-fit 
systemic velocity of 98.54 km s"' corresponds well to that esti- 
mated from molecular lines (Sect.3.1). The SE water emission is 
reasonably well-fitted by two diff'erent models. In the first model, 

01, the vertex of the cone with 20=30° coincides with the centre 
of the methanol ellipse within 25 mas and the velocity of masers 
is 77kms"' (Tab. |5] Fig.|6]l. We argue that this fit is consistent 
with the outflow scenario. 

As we failed to identify any mid-infrared (MIR) counterpart 
to the centre of the methanol ellipse as a powering source of 
the methanol structure and perpendicular to the water maser out- 
flow, we propose an alternative model. This second model, 02, 
assumes that the position of the cone vertex is just between the 



red-shifted and blue-shifted voids of SE water masers. The best- 
fit velocity of the outflow is ~8kms \ while the projection of 
the axis outflow onto the plane of the sky is roughly parallel 
to the main axis of the methanol structure (Tab. |5] Fig. |6). In 
star forming regions, different centres of activity are reported on 
scale sizes of a few to several hundred AUs (e.g., Torrelles et 
al. l200ni20m i. The proper motion studies of G3 1.58 1+00.077 
can clearly provide conclusive evidence for or against this hy- 
pothesis. 

In the case of methanol masers in G33. 641 -00.228, 
we present the two kinematical models for a powering 
source located at the centre of the best-fit ellipse (ARA, 
ADec)=(59 mas, 17 mas) and that located at the MIR position 
(-495 mas, 380 mas). In the first model, Dl, the infall velocity of 
about 1 kms"' is slightly higher than the rotation velocity (Tab. 
HJ. The structure of water masers is fitted by a narrow ouflow of 
20=18° and velocity of 15 km s ' (Tab. |5] Fig. |6]l. It is consistent 
with a jet-like outflow of water masers roughly perpendicular to 
the main axis of methanol structure. In the second model, D2, 
the visible methanol structure is only a small part of the disc ra- 
dius of ~1900AU and the rotation velocity of 3.7 km s"' (Tab. 
nil, which implies that the enclosed mass is about 3OM0. In this 
case, the water masers may trace the wind gas from a disc. 

We note that the kinematic models are unable to accurately 
reproduce the observed structures of the methanol and water 
masers in both targets, and that only measurements of the proper 
motions of maser spots of the two species can unambiguously 
constrain the proposed scenarios. 

5. Conclusions 

We have carried out high-angular resolution studies of the 
22.2 GHz of water maser line towards two methanol maser 
sources G3 1.58 1+00.077 and G33.641 -00.228 using the EVN. 
Although their morphologies did not differ significantly from 
the previous VLA results, the astrometry at the mas level and 
the properties of the maser clusters could be estimated owing to 
the high-angular and spectral resolution. In total, we detected 
eight and six water maser clouds towards G3 1.58 1+00. 077 
and G33. 641 -00.228, respectively. In the first target, the wa- 
ter maser components are associated with different centres of 
star-forming activity, and the components associated with the 
methanol ring-like structure possibly trace the outflow. In the 
source G33. 641 -00.228, southern water masers possibly trace a 
wind from a disc. The kinematic models containing ring-like or 
arc-like methanol maser structures are able to trace a circumstel- 
lar disc/torus around a high-mass young stellar object, whereas 
the water maser distribution is orthogonal to the major axis of the 
applied methanol structure and poorly constrained by the present 
data. The present studies show that the two sources are good 
targets for proper motion studies in order to understand more 
clearly the kinematics of gas in the environments of high-mass 
stellar objects. They also encourage us to extend the multi-epoch 
EVN observations for a whole sample. 
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Table 3. Parameters of 22.2 GHz water and 6.7 GHz methanol maser clouds. 



Cloud Sp Sa, Vfl, FWHM Lp„j Vg,ad Tb 
(mJy beam-') (mjy beam-') (kms"') (kms-') (mas) (AU)'' (km s'^ mas"') (km s"' AU'')' (xlO^ K) 



G31.581+00.077 22.2 GHz water masers 



1 


1008 


969 


102.9 


0.54 


0.17 


0.93 


- 


- 


10.4 


2 


233 


221 


101.8 


0.58 


0.07 


0.38 




- 


2.4 


3 
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444 


100.3 


0.61 


0.47 


2.58 


1.3 


0.24 


4.9 


4 


16534 


- 


- 


- 


0.10 


0.55 


21.2 


3.84 


170.8 


4 1 


- 


5302 


100.2 


0.50 


- 






- 


- 


4 2 


- 


14630 


99.8 


0.67 


- 






- 


- 


4 3 


- 


9445 


99.3 


0.50 


- 






- 


- 


5 


2577 


- 


- 


- 


1.17 


6.42 


1.8 


0.33 


26.6 


5 1 


- 


2137 


97.8 


0.93 


- 






- 


- 


5 2 


- 


1140 


97.4 


0.57 


- 






- 


- 


6 


363 


347 


93.6 


0.72 


0.35 


1.92 


2.0 


0.36 


3.8 


7 
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91.7 


0.94 


0.36 


1.98 




- 
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8 
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1.32 




- 


13.8 












6.7 GHz 


methanol masers 






a 
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99.6 


0.39 


5.8 


31.84 


0.09 


0.016 


3.4 


b* 
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1.1 
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98.8 


0.68 


7.4 


40.63 


0.07 


0.013 


6.4 
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98.7 
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15.92 


0.12 


0.022 
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- 


- 
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f 
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8.78 


0.23 
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- 


- 
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i 
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0.9 


k 
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6.7 GHz methanol masers 
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10.93 


0.21 


0.056 


63.9 


k 


12260 


12665 


59.8 


0.22 


6.0 


22.62 


0.06 


0.016 


38.4 


1 


4105 


4362 


59.6 


0.26 


10.9 


41.10 






12.9 


n* 


109 
















0.3 


P* 


176 
















0.6 


r 


28300 


29569 


58.8 


0.30 


1.2 


4.52 






88.6 



• The emission only in two channels. 

*' More than two spots in a cloud, but no Gaussian characteristic of its spectrum. 

' Assuming the near kinematic distances, 5.5 kpc and 3.8 kpc for 031.581+00.077 and G33. 641-00.228, respectively (Sect. 3). 
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